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Abstract 

We show in detail that the fcT-factorization for exclusive processes is gauge-dependent and incon- 
sistent. 



The /cT^-factorization has been widely used to study exclusive B-decays(see references in [H [2]). In 
this factorization, transverse momenta kx of partons are taken into account. The hard part of the 
factorization is extracted from scattering of off-shell partons and it depends on the transverse momenta. 
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^> , Because the scattering is of off-shell partons, it is likely that the hard part, hence the factorization, is 
I gauge-dependent. To our knowledge, the A^j'-factorization for exclusive B-decays has been not studied at 
^ ■ one-loop completely. 

The first one-loop study of the /cT^-factorization is for the case 7* -|- vr — )• 7 in [1] , where the hard part 
I is calculated in Feynman gauge. It has been claimed that the /cr-factorization is gauge-independent[l]. 
■ In |2j it has been pointed out that the /cr-factorization is in fact gauge-dependent because of singular 
contributions from the wave functions to the hard part in a general covariant gauge. A method has 
been suggested in [3] to eliminate such singular contributions. However, such a method is inconsistent as 
pointed out in [4j. To determine the hard part at one-loop, one needs to calculate at one loop the form 
^ ■ factor in the case and the wave function. It should be noted that there are no complete one-loop results 
i of wave functions in the general covariant gauge. Therefore, the gauge invariance of the /cT^-factorization 
has been never checked explicitly at one-loop with this gauge, except the singular contribution found in 

Recently, the case 7* -|- vr — )■ vr has been studied with the /cr-factorization at one- loop in ^ , where it 
has been pointed out that the gauge invariance of the fc^-factorization is proven in [3] . In this note we will 
study this issue at tree-level and beyond tree-level. We show that the /cy-factorization is gauge-variant 
and hence inconsistent. 

1. Because a gluon is exchanged at tree- level, the problem of the gauge-invariance already appears 
in the A^f-factorization for the form factor of the process 7* -|- vr — t- vr. We use the light-cone coordinate 
system, in which a vector is expressed as a'^ = (a"*" , a~ , a_|_) = ((a" -|- a^) / \/2, (a" — a^) / \/2, a^^a^) and 
Oj^ = (a^)^-|-(a^)^. We take a frame in which the initial- and final vr have the momentum P'^ ~ (P"*", 0, 0, 0) 
and ~ (0,-ftr~,0, 0), respectively. P"^ and K~ are large and the square of the momentum transfer is 
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given by = 2P^K . In the /cr-factorization the form factor for large is written as: 

F^{Q) as J dx(fk± dy(fp_L (l){x,k±)(j){y,p_L)H{x,k_i_,y,p±). (1) 

In the above (j){x, k_i_) and (f>{y,p±) are the wave function for the initial- and final vr, respectively. These two 
wave functions can be different and their definitions will be given later. The hard part H is determined 
by replacing the initial- and final vr with a off-shell quark pair respectively. We take the quark pair 
q{ki)q{P — ki) for the initial tt{P) and the pair q{k2)q{K — ^2) for the final tt{K). The momentum ki 
and /c2 are specified as: 

ki^ = {xoP+,0,ki^), k!^ = {0,yoK-,k2±). (2) 

We have here kf = —kf^ 7^ and /c^ = ~^2± 7^ reflecting the fact that the quark pairs are off-shell. 
To determine H one uses the quark pairs to calculate the form factor and wave functions. For the off- 
shell quark pairs in the initial- and final state one uses the spin projection 757" and 7"^ 75 respectively. 
With these projections one picks up the leading-twist contributions. At tree-level the wave functions, 
denoted as 4>^^\ are proportional to (^-functions and gauge-invariant because no gluon is exchanged. It is 
straightforward to obtain in Feynman gauge H at tree-level, denoted as H^^^ as: 

H^"Hx,k±,y,p±) = ^ . (3) 

xyQ^ + {p±-ki_Y 

In deriving this result one has used the power counting: xQ ~ yQ ^ p± ^ k± ^ 6 and only the leading 
terms in 5 has been taken into account. 

Because it is derived in Feynman gauge with off-shell partons, H can be different in different gauges. 
Supposing we work in the general covariant gauge, in which the gluon propagator reads: 



7^ -|- ie 



+ l£ 



(4) 



with the gauge parameter a, we obtain H as: 



xyQ^ + (p± -k^f\ 2 xyQ^ + {p^_ - ki_f J 

It is clear that H is gauge-dependent at tree-level. The terms with a are at the same order of H^^^ 
determined in Feynman gauge. They are not suppressed by power of 6. Therefore, the gauge-dependent 
term can not be neglected. Similar results for exclusive i?-meson decays are also obtained [6]. Because the 
transverse momenta appear in numerators of the gauge-dependent terms, one may argue that these terms 
may be factorized with higher-twist operators other than the leading-twist operator used to defined (/>'s. 
If one can do so, these terms are still gauge-dependent and can not be neglected with the power counting, 
in comparison with the term factorized with (/)'s. This can be illustrated with the term in the numerator 
which is linearly in k± in Eq.(5). The contribution can be factorized with a wave function defined with 
the matrix element of the initial vr (0|g(0)7^7+9j^g(y)|7r(P)) with y+ = 0. In this case one may need to 
analyze the contribution with the incoming off-shell quark pair combined a off-shell gluon. Adding this 
it may result in that the derivative becomes the covariant derivative Dj_. But it is here not important 
how this term is factorized. The important is that the form factor calculated with off-shell partons has 
a gauge-dependent and nonzero contribution. In Eq.(5) we have factorized it with the wave function. 
Regardless how it is factorized, this gauge-dependent and nonzero contribution can not be eliminated by 
factorization with different operators. Therefore, the /c-r-factorization at tree-level is already inconsistent. 
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An interesting fact should be noted when one studies the /cT^-factorization beyond tree-level in Feyn- 
man gauge. In this gauge, the form factor and the wave functions will not have I.R. or collinear diver- 
gences, because they are regularized by the off-shellness of partons, i.e., by lnA;^_|_ and lnA;2_|_ here. Since 
everything is finite, the factorization is a trivial task. 
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Figure 1: The one- loop diagrams of the wave function. The double line represents the gauge link in Eq.dT]). 

2. To determine the hard part at one-loop, one needs to calculate the form factor and wave functions 
with the quark pairs at one- loop. With these one- loop results and by expanding the right-hand side of 
Eq.(l) in Ug one can obtain H^^\ H^^^ will receive a contribution proportional to 



//(o)^ = dx J d^k±H^°\x,k±,yo,k2±)(l)^^\x,kj 



(6) 



where (j)^^^ is the one-loop contribution to (p calculated with the off-shell quark pair. The wave function 
</> of the incoming vr(P) is defined as: 



{x,kT,C,f^) 



ixP+ z- -iz^-kr 

(27r)3 



(0|g(0)4(oo, 0)7+75L«(oo, z)q{z)\Ti{P)) 



5 + =0 



(7) 



where L„(oo, z) is a gauge link starting from the space-time point z to oo along the vector u. The vector 
u is defined as = (u^, u^, 0, 0). In the above definition the limit ii+ — )• should be taken, i.e., any 
term proportional should be neglected. The limit is taken after all loop-integrations. Besides x, 
kx and the renormalization scale ^ the wave function depends on the vector u through the parameter 
(^^ = 2u^ {P^)"^ /u^ ~ {2u • P)"^ /v?. This results in that the hard part H will also depend on C,. This 
C-dependence is very useful for resummation of double log's. 

Similarly, one can defined the wave function of the outgoing ■n{K) with the gauge link along the 
direction u'^ = (n'"*", 0, 0) in the limit u'~ — t- 0. The wave function of 'k[K) and hence the hard 



Therefore, the two wave 
Q' . Our discussion in the 



part will also depend on C, defined as C'^ = 2u'^{K^Y /u'^ » (2n' • KY /u'"^ 
functions are in general different. In practice one may take the choice C = 
following will not be affected if one makes or does not make the choice. 

The wave function has been studied with an on-shell quark pair at one-loop in The obtained 
results are gauge-invariant. But, the wave function in the /cj^-factorization is calculated with an off-shell 
quark pair, it is not gauge-invariant. Hence, it is possible that the ^-dependence is also gauge-dependent. 
This in turn gives to the hard part H a gauge-dependent ^-dependence, since the form factor does not 
depend on C,. If this is the case, it clearly indicates the gauge-variance of the fer-factorization. In the 
general covariant gauge, because the gauge-dependent term in Eq.(4) is proportional to q^q^ , it is very 
simple to show that this term does not give at one-loop the ^-dependence extra than that in Feynman 
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gauge. However, it is unclear if there exists a gauge-dependent ^-dependence beyond one-loop. The 
situation changes if we work with an axial gauge, there is an extra ("-dependence at one-loop. It is even 
worse that the extra ^-dependence is I.R. divergent. 





Figure 2: The one-loop diagrams of the wave function. The double line represents the gauge link in 



Eq.©. 



We take the axial gauge to examine this. The gauge is fixed hy v ■ G = with the vector 
{v~^ ,v~ ,0,0). and v~ are arbitrary. In the gauge the gluon propagator is given by: 



-|- ie 



^f,u 2 'i_ _^^2 



{v ■ qf 



(8) 



At one-loop, the wave function calculated with the off-shell quark pair is linearly related to the gluon 
propagator because only one gluon is exchanged. Therefore, the possible extra (-dependence can only 
come from the second term in the above. The possible contributions to this are represented by diagrams 
given in Fig.l and Fig. 2. We denote the one-loop contribution from the second term to the wave function 



as 



The calculations are simple, e.g., the contribution from Fig.lc: 



(1) 



-igl5{x — XQ)5^{kT — ^i±)C'f 



Ic 



u 



vTi (7 757+757 • (^1 - 9)7 • q) 



(27r)4 {u-q- ie){q2 + ie){{ki - qf + ie){v ■ q) ' 



(9) 



where we used for the external lines of the off-shell quark pair the spin projection discussed before. 
For the gluon propagator in Fig.lc we only used a part in the second term in Eq.(8). Only this part 
deliveries the extra (-dependence. A caution should be taken for the denominator 1/v ■ q in the above 
when one integrates over q~ . One needs to supply a prescription for the denominator. This can give 
in the integration an additional pole than the pole from the first term. But the contribution from the 
additional pole does not depend on ( in the limit ( — t- 00. In principle there should be a transverse gauge 
link at = 00 in Eq.(7) to make the definition completely gauge invariant. In our case as long as we 
keep and f ^ nonzero, the transverse gauge link will not introduce any contribution. 

It is straightforward to perform the loop integral in Eq.(9) and other loop integrals in Fig.l. In the 
calculation we will meet possible I.R. divergences. We introduce a small mass A for gluons to regularize 
the I.R. divergences. We have from Fig.l: 



C^#(x,A:t,C,m) 
C^cpW^x,kT,C,f^) 



la 



Ic 



Id 



C-^<P^^Hx, kr, C, /i) = <P^''Kx, kT, C, ;u)^ In 4, 



-0W(x,A:t,C,m) 
-</>W(x,A:t,C,m) 



vr 
vr 



, Cx^ 1 , Cx^ 



1 C^(l-x)^ , ^. cHi-x? 

kl +2 //2 



(10) 
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For the contributions from Fig. 2 it is easy to show that the sum does not depend on 



2a 2b 



2c. 



0. 



(11) 



Therefore we have in the gauge v ■ G = Q for the wave function and the hard part H at one-loop: 
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C-^(t>{x,kT,C,f^) = [V'. 
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H{x,kT,C), (12) 



where we use the notation V (p V <Si H to denote the contributions from Feynman gauge, i.e., the 
contributions from the first term in Eq.(8). The difference of signs in the two evolutions reflects the fact 
that the form factor does not depend on From the above results, it is clear that the hard part H, 
hence the /cr-factorization, is gauge-dependent, because it is different in different gauges. In the axial 
gauge, the extra (^-dependence is I.R. divergent because of the contributions from Fig. la and lb. Since 
the hard part receives the C-dependence only from the wave function, the hard part in fact contains a 
C-dependent I.R. divergence. This leads to the conclusion that the factorization is violated in this gauge. 

3. In the general covariant gauge the hard part will receive a soft divergence called light-cone diver- 
gence, as shown in [2]. In [3j a method to eliminate this divergence is suggested. But, this method is 
inconsistent as pointed out in [3]. Here, we explain the inconsistence in detail. 




Figure 3: The contours for the integration of q . 

The gauge-dependent contributions at one-loop are proportional to a in Eq.(4). We denote these 
contributions as (pa- We take Fig.2b as an example. The gauge-dependent part is[2l[3]: 



(j)aix,kT) 



Siaus dq 2{kt - q^)q - • + ?! ^^^^ 



2b 7_oo 27r [{ki - qf + i£\{q^ + ief 



In the above q is the momentum carried by the gluon in Fig. 2b. The components q'^ and q'^ are fixed as 
q'^ = k'^ — xP~^ and gj_ = k\_\_ — with kf = xqP~^. The integral can simply be calculated by taking 
a closed contour in the upper-half- or in the lower-half complex plan of q~, as showing in Fig. 3, where 
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we take the contours consisting of a straight hne along the real axis and a semi circle with the radius R. 
The limit — >• oo should be taken corresponding to that the integration over q~ is from — oo to oo. From 
positions of g~-poles of the integrand the integral only becomes nonzero in the region < < X()P~^, 
because in this region, the pole from the l/[(A;i — g)^ + ie] is in the upper-half plan and the double pole 
from + ie)'^ is in the lower-half plan. For the region < all poles are above the real axis, while 

all pole are below the real axis for q~^ > xqP'^. Taking any one of the two contours the integral can be 
performed. The result does not depend on the choice of contours. 

It is found in [2] that the result is singular. To clearly see this, we take a test function t{x,k±) to 
calculate the convolution t0(f)^'^\ By taking t{x, k±) = H^^\x, k±,yQ, k2±) one obtains the convolution in 
Eq.(6). The singularity can be regularized by a small gluon mass A^, or with dimensional regularization 
as showing in {§]. We first discuss the case with the dimensional regularization, in which the transverse 
space is 2 — ex dimensional with — 0. We have then the singular contribution of (pa\2b and the 
convolution: 



26 



26 



4aa, klJ{q+)e{k+ - q+) 
' kt{ql + q+klJkt)^ 



7rP+ 



+ finite terms. 



t{xo,ki±] 



— ) + (finite terms). 



(14) 



Prom Eq.(13) it is clear that the singularity comes from the momentum region where q has the scaling 
patten q^ ~ ((Jg, 1, (5o, (5o) with 5o ^ 1- The singularity comes from the first term in Eq.(13). We call this 
as light-cone singularity. Adding all contributions at one-loop the singularity is not canceled. It results 
in that the hard part H at one-loop is divergent in the general covariant gauge. It should be noted that 
the form factor does not contain such singular ities [21 19]. 

To deal the singularity a method is suggested in [3]. In the method one keeps R large but finite in 
the calculation of the wave function and the limit 72 — )■ cx) is only taken after that the integrals in the 
convolution are performed. Then one obtains for the wave function not only the contribution from the 
pole but also the contribution from the semi-circle. As suggested in [3| for the contribution from Fig. 2b, 
one takes the contour for < in the upper-half plan and the contour for g+ < in the lower half plan, 
one then has: 



4aa, kl^e{q+)9{k+ - g+) 4aa, 



26 



7r2 ktiq'i + q+klJkty 

TR{e,x,k^) 



e{q+) / d9 + e{-q+) 



J^R{9,x,k±) + 



2{kt -q+){Re 



ie\2 



{2q+Re^<^ - q\Y[2{q+ - kt)Re 



ie 



\k 



i± 



(15) 



where J-r is the first term in Eq.(12) with q~ = Re^^ . There are finite terms denoted as • • •. The terms 
in the second line are from semi circles. If we take the limit i? — )• oo in the above as required in Eq.(13), 
we obtain the result in Eq.(14). This also results in that 4'a\2b is nonzero only in the region 0<q+ <kt 
or < X < xq. As suggested in [3] one should keep R finite here and calculate the convolution first. The 
limit is taken after the integrations in the convolution. In this case one has: 



t< 



26 



^t(xo,/ci^)(^-) + 



- o n I hm 



de 



dq+ + 



de 



dq^ 



t{x,k±)J^R{9,x,k±), (16) 
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Working out the singular part of the integrals in the second line and taking the limit i? — )■ oo one has: 



2b 



8aas,, , . / 1 
— -t Xo,fei_L — 



de 



d9 



t(xo,ki±) ( — ) + (finite terms). (17) 



From the above the convolution calculated with the method from [3j is free from the singularity. Two 
observations can be made at the first look from the above result. Because the limit ii —)• oo is taken after 
the integrals in the convolution, it implies that one introduces a cut-off for q~ with \q^\ < R. It is unclear 
how to implement the cut-off in the definition in Eq.d?]). The finite R results in that the contribution 
is not zero with < corresponding to x > xq. Another observation is that the result depends on 
contours. Different choices of contours give different results. 

For the region with < or q'^ < one always has two possibilities of contours, respectively. In 
the above the contour for < is in the upper-half plan. If we take the contour in the lower-half plan 
with the finite R for < q^ and > q^, we then have: 



iaas k'l^9{q^ 



')e{kt - q+) 4aas 



2b vr2 k+iql+q+kfjkty^ 



de + ei-q^ 



dO 



FR{e,x,ki_) + ■■■,{!%) 



where q^ = q\/{2R). Here the contribution from the double pole starts at = q^ instead of q'^ = 0, 
because the double pole with small enough q'^ can be in the outside of the contour with R. In the above 
the first 0-integral is from 6 = —tt to 9 = because the contour here is in the lower-half plan. Comparing 
with Eq.(15), one realizes that 4>a\2b is different with different contours. This in turn gives different 
results of the convolution. 

Analyzing all possibilities one can have three different results for the following choices of contours in 
the two regions of q~^: (I). All contours are in the upper-half plan or in the lower-half plan. (II). The 
contour for < is in the upper-half plane and the contour for g"'" < is in the lower-half plan. (III). 
The contour for < g"^ is in the lower-half plane and the contour for < is in the upper-half plan. 
The choice (II) corresponding to Eq.(15,16). The different results for these 3 choices are summarized as: 



2b 



2b 



2b 



1 



-t(xo, A;i_|_) — M- (finite terms), for (I), 

8aas 



ttP+ 



t{xo,kij_) - ^^t{xo,ki±_) + (finite terms), for (II), 



t(xo,A:i_L) — + 



1 



8aa, 



t(xQ,ki±) — + (finite terms), for (III). 



(19) 



This is the contour dependence pointed out in [4J. The cancelation of the singularity depends on the 
choice of contours with the method in [3j . 

One may use a small gluon mass A to regularize the singularity. The mass in the general covariant 
gauge is introduced as in [10] : 



+ ia 



q 



2 - A2 ie (g2 - (1 - a)A2 + ie){q^ - A^ + ie) 



(20) 



Comparing with Eq.(4), the double pole in the g~-plan splits into two single poles. It can be shown with 
the method in ^ that the singularity is canceled by the contributions from semi circles. It is interesting 
to note that the cancelation is contour- independent, because an additional contribution from the cases 
where one of the two single poles can be in the outside of a closed contour with the finite R. However, 
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this only works at one- loop, because one has at one- loop only diagrams similar to QED. Beyond one-loop 
level, one can not take a finite gluon mass here for QCD. 

A natural question with the suggest method in [3j is then which answer is correct? When one takes 
the dimensional regularization, the result is contour-dependent. When one uses a finite gluon mass for 
the regularization, it is definitely not correct beyond one-loop level. From the above analysis different 
results about the singularity are due to keeping R finite in Eq.(15) instead of — )• oo. In fact, none of 
them is correct in the sense that R can not be kept finite in the wave function. A finite R is in conflict 
with translational covariance. Below we discuss the conflict in detail. 

We start with the well-known fact that the wave function becomes zero with x > 1. This is derived 
by sandwiching a complete set of states into Eq.© and using the translational covariance. It should be 
noted that xq in Eq.([2]) is arbitrary in the region < xq < 1- For nonzero ki± we can take xq = 1 for 
the off-shell quark pair q{ki)q(P — ki) as an extreme case without any problem, because the quark and 
the antiquark with xq = 1 still have nonzero momenta. Then, with a large but still finite R, the wave 
function is not zero with x > 1. This can also be seen from the contributions in Eq.(15). Therefore, a 
finite R is in conflict with translational covariance in the case with xq = 1. In general, as showing in 
[1], the contributions to the wave function from a class of diagrams including Fig. 2b, where gluons are 
only exchanged between the quark field q{z) and L|j(oo, 0) in Eq.(j7]) or are exchanged to form self-energy 
corrections to q{z) and L]^(oo,0), are proportional to 



2tt (2ttV 



■(0| (4(oo,0) - l) 7+75g(^)k(fci))L+=o (21) 



besides some trivial factors. The leading order result of the above expression is given by Fig. 2b by 
deleting the antiquark line. Sandwiching a complete set of states and using the translational covariance, 
it is easy to show that the above quantity must be zero with x > xq. Here xq is arbitrary in the region 
< Xq < 1. Therefore, the contribution from Fig. 2b must be zero with x > xq. A nonzero contribution 
from Fig. 2b with x > xq is unphysical. Taking \q^\ < R instead of \q~\ < oo in Eq.(13,15), it results in 
that the contribution is not zero with x > xq. One also notes from Eq.(15,18) that with the finite R the 
wave function is not zero for q~^ > kf oi x < 0. It implies that the quark entering hard scattering is with 
a negative energy. This is inconsistent not only with the translational covariance but also with physical 
picture. 

From the above discussion, it is clear that the method of [3] by keeping R finite in the wave function 
is essentially to include the unphysical contributions. Convoluting the wave function containing these 
unphysical contributions with a test function, the discussed light-cone singularities may be canceled. 
The cancelation depends on contours and also on how the singularities are regularized. As shown in the 
above, the existence of these unphysical contributions is in conflict with general principles, i.e., with the 
translational covariance. Because of this and the contour dependence the method is not consistent. We 
emphasize here that the key problem of the method is the inconsistence with the translational covariance 
for the defined wave function in Eq.(7). This in fact can be easily seen by noting the following fact: 
The finite cut-off R for q~ implies that one takes the corresponding space-time coordinate x"*" as one- 
dimensional lattice with the lattice spacing 27r/R. The lattice certainly has no symmetry of translational 
covariance. The conclusion here is, as made in [2j,i4j, that the hard part will contain the divergent part in 
the general covariant gauge. Hence the fcr-factorization is gauge-dependent and violated in this gauge. 

To summarize: With the results presented in this note, one can conclude that the fcr-factorization is 
gauge-dependent and violated in the gauges studied here. Since some sources of the gauge dependence 
studied here are from wave functions and absent in scattering amplitudes of off-shell partons, the fcy- 
factorization is gauge-variant for any exclusive process. Our results can be generalized to the case of 
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-B-meson decays. Therefore, the /cj^-factorization for exclusive i?-meson decays is also gauge-dependent. 
Note added: 

After we completed the present work, an analysis of three-parton contribution to the vr-form factor 
in the A;r-factorization appeared in pjj. The contribution is at twist-three, where the spin projection for 
the initial- and final off-shell quark pair are cr~"75 and cr"'"'^ 75, respectively. In [11] the fc^-factorization 
has also been studied in the general covariant gauge. Because the different spin projections, the results 
obtained in [TT] has no implication for our result in Eq.(5) and our discussion there. But the results with 
the general covariant gauge in [ilj clearly show that the factorization is gauge- variant. 

From the results in [11] the gauge dependent part in the three-parton contribution related to the 
initial vr can be written in the form of the following convolution: 

aH^, (0|g-(y)757+7''^'^(0)g(0)|7r(P)), (22) 

where a is the gauge parameter in Eq.(4) and H^^'^ is perturbative coefficient function. The field q 
is with y'^ = (0,y~ ,y±). It should be noted y± ^ because the quark carries a nonzero transverse 
momentum into hard scattering in the fcr-factorization. For y± = the above matrix-element is of 
twist-3 in collinear factorization. It has been argued in [11] that the above expression is zero because 
of the equation of motion. But it is not true. There are two possibilities to have the above quantity to 
be zero with equation of motion. If one has H^'^ oc g^^ , it is clear that the quantity is zero. In general, 
the perturbative coefficient functions related to the derivative and to the gluon field in the covariant 
derivative D'^{x) = d'^ + igsC^ (x) can be different. In [TT] one has shown that the contribution related to 
the 5"*" can be neglected and the attachment of a -|— component of gluon to the hard scattering diagrams 
gives no contribution. From these results one has = 0. Therefore, H^'^ is not proportional to 

gf^'^ . Another possibility is that the matrix element is proportional to g^'^ . But again this is not correct. 
To show this, one can make a Fourier transformation of (0|g(y)757''7'^-D^(0)g(0)|7r(i-')) and study its 
decomposition. The decomposition is complicated. The ffi'st few terms are: 

= fo{x, k^) {PPg^^-^ - gP^P^) + h{x, k^)gP^^P-^ + f2{x, k^) {gP^k^^ - g^^^k"^) + h{x, k^)gP>^k\^ 
+U{x, k^) {Ppe^kl - P^kP^kl) + • • • , (23) 

where • • • denote the structures which are zero when p = + and those may depend on n^^ = (0, 1,0,0). 
In general the functions /o,i,2,- - are not zero. From the reason of dimension they may be at the same 
level of importance. This can also be seen the constraint of /'s obtained by contracting g^j^y from both 
sides. Taking p = +, one clearly sees that the matrix element is not proportional to g^"^ . Therefore, the 
above quantity is not zero. This shows that the three-parton contribution in the A;y-factorization for the 
TT-form factor in [IT] is gauge-dependent. 
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